Casimir effect between dissimilar materials: a test for the proximity theorem 
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We calculate exactly the Casimir force between a spherical particle and a plane, both with arbi- 
trary dielectric properties, in the non-retarded limit. Using a Spectral Representation formalism, 
we show that the Casimir force of a sphere made of a material A and a plane made of a material B, 
differ from the case when the sphere is made of B, and the plane is made of A. The differences in 
energy and force show the importance of the geometry, and make evident the necessity of realistic 
descriptions of the sphere-plane system beyond the Proximity Theorem approximation. 
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I. INTRODUCTION 

In recent years, accurate experiments IE S IS 
have been performed to measure the Casimir force which 
is one of the macroscopic manifestations of the fluctu- 
ations of the quantum vacuum ,6] . Most of these ex- 
periments H, IE ^ IS have been done using a spherical 
surface and a plane, instead of two parallel planes, as 
originally proposed Casimir 6] . The interpretation of the 
Casimir force in these experiments is based on the Prox- 
imity Theorem which was developed by Derjaguin and 
collaborators to estimate the Casimir force between 
two curved surfaces of radii i?i and i?2- The Proximity 
Theorem assumes that the force on a small area of one 
curved surface is due to locally "flat" portions on the 
other curved surface, such that, the force per unit area is 

where V{z) is the Casimir energy per unit area between 
parallel planes separated by a distance z. In the limit, 
when i?i = R, and i?2 ^ oo, the problem reduces to 
the case of a sphere of radius R and a flat plane, that 
yields to J-{z) — 2-kRV{z). The force obtained using the 
Proximity Theorem is a monotonic function of z, where 
at "large" distances J-{z) (x z~^, while at short distances 
T{z) oc This theorem is supposed to hold when z <C 
i?i,i?2 7 however, the validity of the Proximity Theorem 
approach has not been proved yet. 

Casimir showed that the energy lA (z), between two par- 
allel perfect conductor planes can be found through the 
change of the zero-point energy of the classical electro- 
magnetic field IS, like 

^iz) = ^^[^^^{z)-UJ^{z^ Co)], (1) 
i 

where LiJi{z) are the proper modes that satisfy the bound- 
ary conditions of the electromagnetic field at the planes 
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which are separated a distance z. For real materials, Lif- 
shitz obtained a formula to calculate the Casimir force 
between parallel dielectric planes S]. The force per unit 
area, according to the Lifshitz formula, is 

with ^" = [rf r2 exp(2iA:z)]~^, where r" is the reflection 
amplitude coefficient of the plane j for the electromag- 
netic field with a-polarization. Here, j = 1, or 2, a = s 
or p, q = to / c, Lo is the frequency and c is the speed of 
light. Here, q — (Q, k) is the vacuum wavevector with 
projections parallel to the surface Q, and normal to the 
surface k, k — k + iO+, and g is the photon occupation 
number of the state k. 

The Lifshitz formula depends only on the reflection co- 
efficients of the planes and the separation between them. 
This reflection coefficients can be calculated for arbitrary 
material using the surface impedance formalism for semi- 
infinite slabs @, 0, 0| and for finite slabs using, for 
example, the transfer matrix formalism [T^ . In addi- 
tion, one observes that the Lifshitz formula is symmetric 
under the change of index j = 1 — > 2, it means that 
the force is indistinguishable under the change of planes 
which is natural given the symmetry of the system. 

The force between a sphere and a plane made of ar- 
bitrary dielectric materials is commonly calculated us- 
ing the Proximity Theorem approximation, where the 
Casimir energy V{z) is found using the Lifshitz formula. 
As a consequence, the force between a sphere and a plane 
is indistinguishable if the sphere is made of a material A 
and the plane is made of B, or if the sphere is made of B 
and the plane is made of A. Furthermore, the Proximity 
Theorem assumes that the proper electromagnetic modes 
of the sphere and a plane are the same as those for two 
parallel planes. Such that, the Proximity Theorem ap- 
proach neglects the pressure radiation coming from the 
quantum fluctuations on the opposite side of the sphere. 
In conclusion, in applying the Proximity Theorem one is 
assuming that the geometry of the system is quite irrel- 
evant to flnd the Casimir force. 

In this paper, we show that the geometry of the sys- 
tem is important. By calculating the zero-point energy of 
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the interacting surface plasmons of a sphere-plane con- 
figuration, we obtain that the force between dissimilar 
materials differ if the sphere is made of a material A and 
the plane is made of B, and vice versa. We find a gen- 
eral formula to calculate the energy and Casimir force 
differences between arbitrary materials. We study the 
specific case of aluminum and gold. We restrict ourself 
to the non-retarded limit, this means, to the case of small 
particles at small distances. 



infinite, we obtain the interaction energy of the system 
per unit area, like 
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(6) 



which depends on the energy of the surface plasmon of a 
half space and the distance between planes. 



II. ZERO-POINT ENERGY OF THE 
INTERACTING SURFACE PLASMONS 

A. Parallel planes 

In 1968, van Kampen and collaborators |Q| showed 
that the Lifshitz formula in the nonretarded limit is ob- 
tained from the zero-point energy of the interacting sur- 
face plasmons of the planes. Later, Gerlach ,15] did an 
extension showing that also in the retarded limit the 
Lifshitz formula is obtained from the zero-point energy 
of the interacting surface plasmons. A comprehensive 
derivation of the Lifshitz formula using the surface modes 
is also found in Ref. 

In the case of two parallel planes, in the absence of 
retardation, the proper electromagnetic modes satisfy the 
following expression 



e(w) + 1 



eiiv) - 1 



exp[2fcz] -1 = 



(2) 



for any value of ^ fc ^ 00. The roots uji{k) are identified 
as the oscillation frequencies of the surface plasmons. To 
find the proper modes, it is necessary to choose a model 
for the dielectric function of the planes. To illustrate the 
procedure for metallic planes, let us employ the Drude 
model 



e{Lu) = 1 - 



Lu{uj + i/r) ' 



(3) 



with LOp the plasma frequency and r the relaxation time. 
Substituting Eq. in Eq. |j2Jl, and taking into account 
that for most materials 1 S> (rwp)"^, we find two proper 
modes that, at large distances (kd :s> 1), are given 



UJ± 



^ (^1 ± i exp[-fcd] - i exp[-2fcd] ± • • , (4) 



and in the limit kd — > 00, we recover uj± = uj^j^f^ which 
is the frequency of the surface plasmon of each half space. 
The zero-point energy per unit area will be given by 



-1 T- /"CXJ 



(5) 



Now, taking the difference of the zero-point energy when 
the planes are at a distance z, and when they are at 



B. Sphere and plane: dipolar approximation 

Now, let us consider the case of a sphere of radius a and 
dielectric function 6^(0;), such that, it is at a minimum 
distance z from a plane of dielectric function £^(0-"), as 
shown in Fig. 1. The quantum vacuum fluctuations will 
induce a charge distribution on the sphere which also 
induces a charge distribution in the plane. Then, the 
induced /m-th multipolar moment on the sphere is given 
by [13 

QUu^) = + (^)] , (7) 

where VJ™'^(i-^) is the field associated to the quantum vac- 
uum fluctuations at the zero-point energy, Vf^iyi) is the 
induced field due to the presence of the plane, and ai„x(LS) 
is the /m-th polarizability of the sphere. To show the 
procedure in detail, first we work in the dipolar approxi- 
mation, i.e., when / = 1. In section III, we will calculate 
the proper modes exactly, taking into account all the 
high-multipolar charge distributions. 



FIG. 1: Schematic model of the sphere-plane system. 

Using the method of images, the a relation between the 
dipole moment on the the sphere Ps{uj), and the induced 
dipole moment on the plane 'Pp{uj), is 



Pp(w) 



1 - ep(^) j 
1 -I- ep(w) 



•Ps(^)- 



(8) 



Here, M = (—1,-1,1) is a diagonal matrix whose ele- 
ments depend on the choice of the coordinate system. 
From Eq. {Tj), one finds that the total induced dipole 
moment on the sphere is 



(9) 



where the interaction between dipoles is coupled by the 
dipole-dipole interaction tensor 

Here, I is the identity matrix, and r is the vector be- 
tween the centers of the sphere and the dipole-charge 
distribution on the plane. From Fig. 1, we find that f = 
(0,0,2(z-|-a)), such that, M-T = (-l/r^, -l/r^, -2/r3) 
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is also a diagonal matrix. Substituting Eq.® in Eq.®, 
one finds 



1 



a{uj) 



I + fc{uj)M ■ T 



(10) 



where G(cj) is a diagonal matrix with 

/eH = [l-epH]/[l+6pH]. 

Multiplying Eq. (|1U|I by a^, one finds that each term of 
G(w) is dimensionless and has two parts: the left-hand 
side which is only associated to the material properties of 
the sphere through its polarizability l/a{oj) = a^/a{uj), 
and the right-hand side which is related to the geometri- 
cal properties of the system trough a, and z. The right- 
hand side also depends on the material properties of the 
plane trough the function fd^)- 

The eigenfrequencies of the sphere-substrate system 
must satisfy Eq. (|10|l . and are independent of the ex- 
citing field, in this case E'^^'^[uj). Then, these frequencies 
can be obtained when the determinant of the matrix in 
the left-side of Eq. H1U|I is equal to zero, detG(ti') = 0. 
Thus, the eigenfrequencies are given by 



1 



1 



2/c(w)a3 
^ [2(2 + a)]= 



0. (11) 



Until this point, it is necessary to consider a model for the 
dielectric function of the sphere to find the proper modes. 
Again, we illustrate the procedure for metallic sphere and 
plane, using the Drude model of Eq. (|3J). Moreover, we 
use the fact that the polarizability of the sphere, within 
the dipolar approximation and in the quasi-static limit, 
is given by 



a{uj) 



i es(cj) ^ 1 
es(w) + 2' 



surface plasmon of the plane and we denote them with 
. These modes are blue-shift as the sphere approaches 
the plane, also because of the interaction between the 
sphere and the plane. In this case, when the separation 
z/a increases, the modes uj- go to Wp/\/2 which is the 
surface plasmon of the isolated plane. This behavior of 
the modes is clearly observed in Fig. ^jp, where uj+ and 
LU- are plotted as a function of z/a. This shows that 
the the zero-point energy of the sphere-plane system is 
directly associated to the interacting surface plasmons of 
the sphere and the plane. 



III. EXACT CALCULATION OF THE 
ZERO-POINT ENERGY OF A SPHERE ABOVE 
A PLANE 



In this section we calculate exactly the proper electro- 
magnetic modes of the sphere-plane system, within the 
non-retarded limit, by including all the high-multipolar 
interactions. These proper modes are calculated using a 
Spectral Representation formalism 18, 19], then, we cal- 
culate the zero-point energy using Eq. QJ. The Spectral 
Representation (SR) formalism has the advantage that 
the contributions of the dielectric properties of the sphere 
can be separated from the contributions of its geometri- 
cal properties. The latter allows to perform a systematic 
study of the sphere-plane system that may be very help- 
ful for experimentalists. The details of the SR formalism 
can be found in Refs. El and El here we only explain it 
briefly. 

We now calculate the proper modes of the system in- 
cluding all the high-multipolar charge distributions due 
to quantum vacuum fluctuations. Using the method of 
images, we flnd for the ^m-th multipolar moment on the 
sphere, that 



FIG. 2: Plot of l/a(tj) in solid line and /c(w)[2(l + z/a)]'^ 
as a function of uj/uip for different values of z/a, for a gold 
sphere over a gold plane. 



I'm' 



AttSu/S„ 



{2l + l)ai>m'{L0) 



+ /c(^)4 



Q.rn'H^mrico), 

(12) 

where A^™ is the matrix that couples the interaction be- 
tween sphere and substrate E3j and aim(w) is the Zm-th 
polarizability of the sphere. The expression for the in- 
In Fig.ISi, we plot, l/&{uj) and /c(c^)[2(l + z/a)]- as teraction matrix A^"" is given in the Appendix section. 



a function of uj/u>p for different values of z/a, for a sphere 
and a plane both made of gold. The proper frequencies 
of the system are given when the black-solid line and the 
color-dotted lines intersect each other. We observe that 
there are two different proper modes for each term on 
Eq. Hll|l . giving a total of six modes for each z/a. The 
proper modes at the left-hand side in Fig. |3i correspond 
to the surface plasmon of the sphere that we denote with 
These modes are red-shift as the sphere approaches 
the plane because of the interaction between them. When 
the separation z/a increases, the modes uj+ go to uip/^/S, 
which is the surface plasmon of the isolated sphere. The 
proper modes at the right-hand side correspond to the 



and it shows that depends only on the geometri- 

cal properties of the system. If we have a homogeneous 
sphere, its polarizabilities are independent of the index 
m, and are given by E3 



ai{uj) 



Ijesjoj) - 1] 
l[es{uj) + 1] + 1 



,2/+l 



nio 



nio - u{u!) 



-,21 + 1 



(13) 



where nio = l/{2l + 1), and u(w) = [1 — €s{uj)]~^. Notice 
that also in the latter equation the dielectric properties 
of the sphere are separated from its geometrical proper- 
ties. By substituting the right-hand side of the above 
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equation, we rewrite Eq. H12() as: 



(;^2/+l)l/2 



47r 



^rH, (14) 



where the /m, I'm'-th element of the matrix H(w) is 

(15) 

The proper modes are given when the determinant of the 
quantity between parenthesis in Eq. (|14() is equal to zero, 
det [—^(0-1)1 + ]HI(cj)] = or when the product 



\][-ti{uj) + n,{uj)] =0, 



(16) 



where ns{u}) are the eigenvalues of IHI(w) for each z. 

The Spectral Representation formalism can be apply 
in a very simple way. First, we have to choose a substrate 
to calculate the contrast factor /c(w), with this, we con- 
struct the matrix for a given z, and we diagonalize 
it to find its eigenvalues Us {uj) . We can repeat these steps 
for a set of different values of z. Notice that the eigenval- 
ues of can be found without to do any assumption 
on the dielectric function of the sphere. Once we have the 
eigenvalues as a function of z, finally, we have to consider 
an explicit form of the dielectric function of the sphere, 
such that, we calculate the proper electromagnetic modes 
ujsifjj) trough the relation u(tt's) — ns{uj). 

To illustrate the procedure, we use again the Drude 
model for the dielectric function of the sphere, and a 
plane with a constant dielectric function. In this case, 
the proper modes are given by 



UJsiz) = -t/T + ^-(1/t)2 + UJ^pUsiz), (17) 

with s — (l,m), and where the eigenvalues ns{z), are 
independent of the frequen cy. If we consider that 1 ^ 
{Tujp)~^, then LOs{z) w ns{z). In this case, the zero- 
point energy, according with Eq. is given by 



(18) 



l.m 



where Hup^/nio = hup^/l/ {21 + 1), are the energies as- 
sociated to the multipolar surface plasmons resonances 
of the isolated sphere. Here, , I = 1,2, ...,L, where 
L is the largest order in the multipolar expansion. For 
I — 1 we have that hwp^JniQ = hujp/\/3, that is the 
surface plasmon of the sphere in the dipolar approxi- 
mation. When I ^ 1, we have that the surface plas- 
mons of the sphere due to high-multipolar moments are 
very similar to the surface plasmons of the plane, because 
hwp^TLifj — !■ hhjp/ \/2 when ^ — > 00. In the presence of the 
plane, the proper modes of the isolated sphere are mod- 
ified according with Eq. H15|) . These proper modes are 



red-shift always as the sphere approaches the plane, and 
this shift depends on z/a in a non-monotonic way. As 
z/a^ Q, more and more multipolar interactions must be 
taken into account. Since the proper modes are red-shift 
as z/a — > 0, it might be possible that ns{z) ^ (TWp)"^. 
In such case, we will have to consider also temperature 
effects, because it is possible to excite modes of very low 
frequencies jjfl]. In this work, we restrict our selves to 
distances such that ns{z) ^ (rwp)"^, temperature effects 
will be studied elsewhere. 

It is possible to find the behavior of the energy as a 
function of z/a for any Drude's sphere independently of 
its plasma frequency, like S = £/fiWp. Furthermore, the 
force can be also studied independently of the radius of 
the sphere, since £{z/a), one finds that 



F 



d£{z/a) d£{z/a) d{z/a) 



dz 



d{z/a) dz 



(19) 



and then, we can define also a dimensionless force like 
aF = aF/hup. In summary, we have shown, for a Drude 
sphere with 1 ^ (rwp)~^, that given the dielectric prop- 
erties of the substrate the behavior of the energy and 
force is quite general, showing the potentiality of the 
Spectral Representation formalism. 

In Fig. 13 we show the dimensionless energy, £, and 
force aF, for a sphere and a plane both made of gold. We 
show £ and aF as a function of z/a when all multipolar 
interactions are taken into account, as well as only dipo- 
lar, and up to quadrupolar interactions are considered. 
The largest multipolar moment to achieve convergence 
of the energy at a minimum separation of z/a = 0.1 is 
L = 80. We observe that £ and aF are a non-monotonic 
functions of z /a when we consider multipolar interaction 
greater than dipolar ones. Only within the dipolar ap- 
proximation, £ and Fa behaves like {z/a)~^, and {z/a)~^ 
for any z, respectively. Let us first analyze the sys- 
tem when only up to quadrupolar interactions are taken 
into account. In that case, the energy as well as the 
force, show three different regions: (i) at large distances, 
z > 5a, only dipolar interactions are important, (ii) 
when 5a > z > 2a, we found that dipolar-quadrupolar 
interactions become important and £ and Fa behaves 
like (z/a)"'*, and {z/a)~^; while (iii) at small distances 
z < 2a the quadrupolar-quadrupolar interaction domi- 
nates, and £ and Fa behaves like (z/a)~^, and (z/a)~^, 
respectively 



FIG. 3: (a) Energy £, and (b) force aF as a function of z/a, 
with L = 80 (dotted line), L = 2 (solid line), and L = 1 
(dashed hne). (c) Force difference, \{F^" - F^^)/F^"\, us- 
ing multipolar moments LH = 80 and LW = 2 (solid line), 
LH = 80 and LW = 1 (dashed line), finally LH = 2 and 
LW = 1 (dot-dot-dashed line). 

However, as the sphere approaches the substrate the 
interaction between high-multipolar moments becomes 
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more and more important. The attractive force sud- 
denly increases more than four orders of magnitude, as 
compare with the dipolar and quadrupolar approxima- 
tions. In Fig. 121; we show the difference of the dimcn- 
sionless force, \{F^" ~ F^^)/F^"\, as a function of the 
separation, with LH and LW the highest and lowest- 
multipolar moments taken into account. We observe that 
at distances z > 2a the force can be obtained exactly 
if up to quadrupolar interactions are considered. We 
also found that for z > 7a the interaction between the 
sphere and the substrate can be modeled using the dipo- 
lar approximation only, like in the Casimir and Polder 
model ^j]. However, at separations smaller than z < 2a, 
the quadrupolar approximation also fails, and it is nec- 
essary to include high-multipolar contributions. For ex- 
ample, at z = a the quadrupolar approximation gives 
an error of about 10%, while the dipolar approximation 
gives an error of about 40%. At smaller distances like 
z = a/2, the quadrupolar approximation gives an error 
of about 40%, while in the dipolar approximation the er- 
ror is 70%. At z = 0.1a both approximations give an 
error larger than 90%. To directly compare the dipolar 
and quadrupolar approximations, we also plot in Fig. 13; 
the force difference when LH = 2 and LW = 1. In this 
curve, we clearly observe that even at z = 5a there are 
differences between the dipolar and quadrupolar approx- 
imations of about 5%. 

The large increment of the force at small separations 
due to high-multipolar effects could explain the physical 
origin of, for example, the large deviations observed in 
the deflection of atomic beams by metallic surfaces p^ . 
as well as some instabilities detected in micro and nano 
devices. Furthermore, to compare experimental data of 
some experiments |^ |j, Q with the Proximity Theorem 
approximation, it has been needed to make a significant 
modification to the Casimir force, since they measured a 
larger attractive force than the one predicted. They at- 
tributed the deviations from the Proximity Theorem ap- 
proximation to the roughness of the surface which tend 
to increase the attraction force, essentially as a series of 
the inverse powers of the separation, z. These devia- 
tions might be attribute also to the interactions between 
high-multipolar charge distributions due to quantum vac- 
uum fluctuations. In summary, the measurements indi- 
cate that dispersive forces between a polarizable atom or 
a spherical particle and a planar substrate involves more 
complicated interactions than the simple dipolc model of 
Casimir and Polder or the Proximity Theorem approxi- 
mation. 

In this paper, we propose a way to elucidate experi- 
mentally two important issues related with the Casimir 
force: (i) the relevance of the geometry via the surface 
plasmon interactions between bodies, and (ii) the rele- 
vance of the high-multipolar interactions in the sphere- 
plane model. In the next section, we study the case of 
the casimir force between a sphere and a plane of dissim- 
ilar materials that can help us to clarify these important 



IV. CASIMIR FORCE BETWEEN DISSIMILAR 
MATERIALS 

Using the formalism presented in the previous section, 
we calculate the energy and casimir force for a sphere 
and a plane of dissimilar materials. To illustrate our re- 
sults, we choose gold and aluminum materials described 
by the Drude model of Eq. (PI, with plasma frequen- 
cies hujp = 8.55 eV, and 15.8 eV, and relaxation times 
(rwp)-! = 0.0126, and 0.04, respectively. 

In Fig. 21 we show the exact calculated energy and 
force for an aluminum sphere above a gold plane, f ai/Au 
and ai^Ai/Auj both in solid line, and the same for a gold 
sphere above an aluminum plane, ^au/ai and ai^Au/Ah 
both in dashed line. In both figures, we find that the 
differences are small and cannot be directly observed. In 
Fig-Et, we show the differences of energy and force given 
by 



and 



A£ = 2 



AT ^2 



E 



A/B 



& 



B/A 



e 



A/B + ^B/A 



A/E 



^1 



B/A 



A/E 



B/A 



(20) 



as a function of z/a, with A=A1, and B=Au. Here, we 
observe that for these materials the differences of energy 
and force vary as a function of z/a. At very small dis- 
tances, (z < 0.5a) the difference of energy and force, un- 
der the interchanged of materials between the sphere and 
the plane, is less than 3%. Then, the difference suddenly 
increases with the minimum separation, and at z > a, the 
difference of energy and force have almost reached their 
maximum. For these particular materials, the maximum 
difference in energy and force is about 6%, independently 
of the radius of the sphere. 



FIG. 4: Exact (a) energy and (b) Casimir force as a function of 
zja, for a sphere made of AI over an An plane (solid line), and 
for a sphere made of Au over an Al plane (dashed line), (c) 
Energy (solid line) and force (dashed line) differences between 
a sphere made of Al over an Au plane , and for a sphere made 
of Au over an Al plane. 

In the previous section, we show that for a minimum 
separation larger than 2a, the force can be described with 
quadrupolar interactions between the sphere and plane. 
This means that for z > 2a the energy and force are 
monotonic functions of z. This conclusion can be easily 
deduce from Fig.O;. We also show in section II, that the 
Casimir energy is given by the interaction of the surface 
plasmon of the sphere and plane. In Fig. |2t>, we observe 
that the proper modes of the system for z > a are ap- 
proximately given by fihjp/y/2 and hujp/^/^, when L = 1. 
These yield to conclude that for z > a, the energy and 
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force differences are quiet independent of z, and are like 



73 ^2 



V3 





tile's \ 






V V3 ' 









= 0.202 



(21) 



For the case of A=A1, and B=Au, we find that AjF ~ 
6.0%, which is the value shown in Fig. ^ when z > 2a. 
Then, as the materials become more dissimilar the dif- 
ference of energy and force become larger, for exam- 
ple, for potassium and aluminum we would find that 
AJ^ ~ 12.4%, while for gold and copper, AJ' ~ 2.4%. 
When z < a the difference decreases as z also does. As we 
see previously, as the sphere approaches the plane more 
and more multipolar interactions become important in 
the calculation of the Casimir force. The energy asso- 
ciate to the interacting high- multi polar surface plasmons 
is proportional to hwp^/l/ (21 + 1), and when I increases, 
hujp^l/ [21 + 1) — > Swp/-\/2, the energy associate to the 
surface plasmon of a semi-infinite plane. Therefore, as 
the sphere approaches the plane, the interaction between 
them resembles the interaction between a plane and a 
set of planes, where each one of they contributes to the 
energy with a factor of (z/a)~", with 3<n^2L+la 
positive integer, where L is the largest order of the in- 
teracting multipolar surface plasmon. In conclusion, as 
the sphere approaches the plane, the interaction between 
them is like the interaction between a plane and a set of 



planes, and not like the interaction between two planes, 
as it is found in the Proximity Theorem approximation. 



V. SUMMARY 

In this work, we show that the Casimir force between 
a sphere and a plane can be calculated through the en- 
ergy obtained from their interacting surface plasmons. 
Our result is in agreement with the work of van Kampcn 
and collaborators 131 , and Gerlach that showed that 
the Casimir force between parallel planes can be obtained 
from the zero-point energy of the interacting surface plas- 
mons of the planes. Then, we showed that the Casimir 
force of a sphere made of a given material A and a plane 
made of a given material B, is different from the case 
when the sphere is made of B, and the plane is made of 
A. We obtained a formula to estimate the force difference 
for Drude-like materials. We found that the difference de- 
pends on (i) the plasma frequency of the materials, and 
(ii) the distance of separation between sphere and plane. 
We also found that as the sphere approaches the plane, 
the force difference becomes smaller and resembles the 
interaction of a plane and a set of planes. 
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APPENDIX 

The tensor that couples the interaction between the sphere and the plane is 0| 
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where Yim(0,ip) are the spherical harmonics. In spherical coordinates r = (2(z + a), 0, 0), then, the spherical harmonics 
reduce to |23l 
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This means that m = m and multipolar moments with different azimuthal charge distribution are not able to couple 
or interact among them. This yields to 
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which is a symmetric matrix. 
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